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Abstract 


Synthetic  aperture  radar  (SAR)  engineers  continue  to  search  for  algo¬ 
rithms  and  techniques  that  reduce  sidelobes  and  improve  the  resolution 
of  a  given  data  set.  The  classical  estimation  of  the  power  spectrum  of  a 
signal  from  a  finite  record  of  data  is  typically  made  using  a  discrete 
Fourier  transform  (DFT).  To  mitigate  errors  due  to  the  data  record  being 
finite,  one  typically  applies  various  windows  (or  apodizations)  before 
transformation.  These  windows  fall  into  two  categories:  raised  cosine- 
based  (e.g.,  Hanning,  Hamming,  Blackman-Harris,  and  Nutall  windows), 
and  noncosine-based  (e.g.,  Dolph-Chebyshev,  Taylor,  Kaiser,  and 
Gaussian  windows).  These  windows  trade  off  resolving  power  for  im¬ 
proved  (reduced)  sidelobes.  Similarly,  SAR  imagery  exhibits  sidelobes 
because  the  images  are  derived  from  a  finite  aperture.  The  same  windows 
are  applied  to  SAR  data  to  reduce  the  sidelobes.  This  report  describes  a 
method  of  allowing  a  raised  cosine-based  window  to  be  adaptively 
changed  at  each  output  data  point  (or  pixel).  The  technique  is  spatially 
variant  apodization  (SVA).  This  spatially  adaptive  window  maintains  the 
resolution  normally  associated  with  rectangular  weighting.  However,  it 
simultaneously  reduces  the  sidelobes  commensurate  with  the  order  of  the 
filter  (i.e.,  the  number  of  cosine  terms  used).  Results  are  shown  on  imag¬ 
ery  from  the  U.S.  Army  Research  Laboratory's  ultra-wideband 
Boom-SAR  system.  It  is  important  to  note  that  SVA  is  not  limited  to  SAR 
applications.  It  is  applied  in  the  image  domain  and  is  applicable  to  all 
systems  that  produce  images. 
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1.  Introduction 

Synthetic  aperture  radar  (SAR)  integrates  energy  over  a  finite  aperture  to 
create  pixel  values  in  the  output  image.  For  a  uniform  (rectangular) 
weighting  function  across  the  aperture,  the  impulse  response  (IPR)  for  a 
point  target  return  in  one  dimension  is  a  sine  function.  The  closest 
sidelobe  is  13.5  dB  below  the  peak,  and  the  sidelobe  envelope  decreases 
6  dB/octave.  Because  the  ultra- wideband  (UWB)  SAR  imagery  can  have  a 
large  dynamic  range,  the  sine  function  sidelobes  from  a  strong  target  can 
easily  obscure  weaker,  nearby  targets.  To  mitigate  errors  due  to  the  data 
record  being  finite,  one  typically  applies  various  windows  (or  apodiza- 
tions)  before  transformation.  These  windows  fall  into  two  categories: 
raised  cosine-based  (e.g.,  Hanning,  Hamming,  Blackman-Harris,  Taylor, 
and  Nutall  windows  [1]),  and  noncosine-based  (e.g.,  Dolph-Chebyshev 
and  Kaiser  windows).  With  the  Hamming  window,  for  example,  the  peak 
sidelobe  is  43  dB  below  the  peak  of  the  response.  This  lower  sidelobe 
comes  at  the  cost  of  a  main  lobe  that  is  1.36  times  wider,  as  compared  to 
uniform,  or  rectangular,  weighting. 

Classical  super-resolution  techniques,  such  as  Capons'  minimum  variance 
method  (MVM),  ESPRIT,  and  MUSIC  (multiple  signal  classification),  are 
based  on  an  assumed  model  and  estimates  of  the  parameters  of  the  model 
given  the  observed  data.  The  solution  of  the  system  of  equations  required 
in  these  techniques  is  computationally  burdensome.  Compared  with 
MVM,  spatially  variant  apodization  (SVA)  does  not  require  inversion  of 
the  covariance  matrix.  Also,  estimating  model  orders  is  not  trivial,  and 
SVA  has  none  to  estimate.  Jung  and  Munson  [2]  have  shown  that  one¬ 
dimensional  (1-D)  SVA  is  a  special  version  of  the  MVM. 

SVA  is  the  technique  that  optimizes  the  IPR  on  a  point-by-point  basis  to 
minimize  sidelobes  in  the  SAR  imagery.  Stankwitz  and  others  [3,4]  intro¬ 
duced  the  1-D  SVA  in  1994.  It  is  an  extension  of  dual  apodization,  which 
generates  two  images  using  two  different  normalized  windows,  and  takes 
the  output  pixel  from  whichever  image  has  the  minimum  absolute  value. 
In  this  report,  we  extend  these  results  to  a  full  2-D  optimization.  We 
describe  several  techniques  that  obtain  high  resolution  and  low  sidelobes. 
They  are  simple  to  compute  and  do  not  assume  a  model.  Instead,  they 
allow  a  sum-of-raised-cosines-based  window  to  be  adaptively  changed  at 
each  output  data  point  (or  pixel)  to  minimize  sidelobe  responses.  This 
adaptive  window  maintains  the  resolution  normally  associated  with 
rectangular  weighting,  yet  simultaneously  reduces  the  sidelobes  com¬ 
mensurate  with  the  order  of  the  window  (i.e.,  the  number  of  cosine  terms 
used). 

Historically,  SVA  has  referred  to  the  special  case  where  the  window  is  a 
single  raised  cosine  and  where  the  cost  function  is  minimum  energy  in 
the  output  pixel  given  a  constrained  region  of  weights  spanning  from 
uniform  weighting  to  a  Hanning  window.  Given  the  derivation  of  the 
SVA,  DeGraaf  [5,6]  introduced  the  adaptive  sidelobe  reduction  (ASR) 
algorithm.  It  uses  the  more  general  sum-of-raised-cosines  window,  and  a 
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cost  function  constraint  to  handle  the  overdetermined  parameter  space. 
That  cost  function  simply  constrained  the  energy  in  the  weights.  A  tech¬ 
nique  to  apply  SVA  multiple  times  to  obtain  super-resolution  has  also 
been  described  by  Stankwitz  and  Kosek  [7], 

None  of  the  previous  papers  addresses  using  the  algorithm  on  ultra- 
wide-bandwidth/ultra-wide-angle  (UWB/UWA)  SAR  data.  The  deriva¬ 
tions  were  based  on  narrow  angle  assumptions  (i.e.,  down-range  and 
cross-range  sidelobes  are  orthogonal  and  lie  on  a  Cartesian  grid).  We 
found,  however,  that  the  technique  can  be  applied  to  UWB/UWA  data 
with  90°  integration  angles.  More  recently,  Rau  and  McClellan  [8]  have 
shown  that  applying  a  wedge-shaped  window  to  the  2-D  Fourier  trans¬ 
form  of  a  UWB  image  is  nearly  equivalent  to  applying  a  sliding  (constant 
integration  angle)  window  to  the  raw  UWB/UWA  SAR  data.  The  result  is 
valid  for  constant  integration  angle,  strip-map  mode  focusing  processes. 
In  other  words,  a  focusing  process  that  produces  an  image  with  a  spa¬ 
tially  invariant  IPR.  The  approximation  used  to  obtain  this  result  is 
accurate  for  the  large  90°  integration  angles  used  in  foliage  penetration 
(FOPEN)  SAR.  Therefore,  optimizing  SVA  to  the  UWA  case  requires  the 
SVA  sample  spacing  to  be  wedge-shaped  (see  eq  (3)).  Fixing  R  for  small 
image  chips  has  been  found  to  produce  excellent  results. 

Figure  1  shows  the  operations  involved  with  image  formation  and  how 
SVA  is  applied.  An  aperture  weighting  is  typically  applied  in  the  raw  data 
domain,  as  shown  in  the  Aperture  window  1  box.  Instead  of  applying  it 
there,  Rau  and  McClellan  showed  that  an  angle-space  window  could  be 
applied  in  the  ( kx ,  ky)  domain,  shown  as  the  Aperture  window  2  box.  On 
the  other  side  of  the  Fourier  transform,  this  angle-space  window  can  be 
applied  using  convolution  instead  of  multiplication  in  the  frequency 
domain.  This  position  is  naturally  suitable  for  pixel-by-pixel  optimization 
because  forward  and  inverse  Fourier  transforms  are  unnecessary  to  test 
different  windows.  Furthermore,  this  position  allows  SVA  to  be  applied  to 
any  image,  regardless  of  how  it  was  made.  Therefore,  it  can  be  used  to 
benefit  more  than  just  the  SAR  community. 
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2.  SVA  Algorithm  Description 

First,  we  must  define  the  sum-of-raised-cosines  window.  Let  N  be  the 
number  of  points  in  an  aperture,  zv  be  the  variable  parameter  per  second 
for  the  window,  and  n  =  -N/2. . .  N/2  -  1.  The  raised  cosine  weighting 
functions  are  given  by 

A(n)  =  1  +  2w1  cos  +  2w2  cos  +  ...  +  2wk  cos  |^^-)  .  (1) 

For  example,  the  family  of  first-order  weighting  functions  ranges  from 
uniform  (w1  =  0)  to  Hamming  (zv1  =  0.426)  to  Hanning  (zv1  =  0.5).  A 
Blackman-Harris  window  is  second-order  with  zv1  =  110/189,  and  zv2 
=  715/7938. 

The  SVA  algorithm  is  derived  by  the  following  argument.  First,  with  the 
raised-cosine  windows  defined  as  above,  the  integral  or  summation 
across  the  aperture  weights  is  independent  of  all  zv.  Therefore,  a  unit  point 
scatterer  (amplitude  =  1),  will  always  focus  to  the  same  amplitude,  re¬ 
gardless  of  all  zv.  This  means  that  the  zv's  can  be  adjusted  to  minimize  the 
amplitude  of  an  output  pixel  because  the  zv's  cannot  affect  the  amplitude 
of  the  desired  point  scatterer.  The  result  is  that  sidelobes  from  other 
targets  (i.e.,  at  positions  other  than  the  desired  output  data  point  or  pixel) 
are  minimized. 

For  example,  if  we  limited  ourselves  to  a  finite  set  of  M  zc-vectors,  SVA 
would  be  equivalent  to  focusing  the  image  M  times,  once  for  each 
re-vector,  and  then  letting  each  pixel  in  the  output  image  equal  the  mini¬ 
mum  (for  the  same  pixel-position)  from  the  M  images  produced  by  the  M 
apodizations.  Such  a  mechanization,  of  course,  would  require  too  much 
computation  to  be  practical. 

To  make  the  approach  practical,  we  first  recall  that  multiplying  by  the 
raised  cosine  window  in  the  (kx,  ky)  domain  (i.e.,  the  2-D  Fourier  trans¬ 
form  of  the  image)  is  equivalent  to  doing  convolution  in  the  image  do¬ 
main.  Taking  the  length-N  the  discrete  Fourier  transform  (DFT)  of  the 
raised  cosine  windowing  function  A(n)  yields  the  Nyquist-sampled  IPR: 

a(m)  =  zvk8m  _k+...  +  w1 8m  +  8ni/  0  + zvxSmA+  ...  +  zvk8nh  k  ,  (2) 

where  8m  k  is  the  Kronecker  delta  function  8m  k 

Thus,  a  first-order  raised  cosine  windowing  function  IPR  contains  only 
three  nonzero  terms;  namely  a(m)  =  zol8m  +  8m  0  +  zv-l8m  v  This 
windowing  function  can  be  implemented  efficiently  using  a  three-point 
convolution  in  the  image  domain  on  complex  Nyquist-sampled  imagery 
that  was  generated  with  a  uniform  window.  If  the  image  is  oversampled 
at,  say,  R  times  the  Nyquist  rate,  then  equation  (2)  becomes 

a(m)  =  zvk8m/ _kR  +  . . .  +  zv ,8^ _R  +  8m  0  + zv18m  R+ ...  +  zvk8nh kR  .  (3) 


1 1,  m  =  k 
|  0,  m^k 
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R  in  equation  (3)  can  be  made  proportional  to  range  to  produce  the 
wedge-shaped  window  described  by  Rau  and  McClellan  [8]. 

Next  we  must  compute  the  optimum  weights  for  each  output  point  (or 
pixel).  The  1-D  SVA  technique  above  must  be  modified  to  apply  to  the 
2-D  SAR  imagery  The  SVA  can  be  applied  in  two  ways:  sequentially  in 
1-D,  first  on  one  axis  and  then  upon  the  other,  or  in  2-D  simultaneously 
The  sequential  approach  has  the  disadvantage  of  not  being  able  to  guar¬ 
antee  that  the  data  provided  to  the  second  application  of  the  SVA  form  an 
unambiguously  Nyquist-sampled  image  because  some  points  have  been 
set  to  zero.  Therefore,  it  may  no  longer  be  a  smooth  or  band-limited 
function.  The  effect  was  used  by  Stankwitz  and  Kosek  [7]  to  obtain  super¬ 
resolution. 

Application  of  the  SVA  on  the  two  dimensions  simultaneously  uses  the 
following  2-D  IPR: 


wmwn 

wn 

wmw 

wm 

1 

Wmwn 

zvn 

wmw 

(4) 


We  will  now  derive  the  optimum  weights.  Let  the  complex  pixel  value  of 
the  original  unweighted  image  be  defined  as 

g(m,  n)  =  gr(m,  n)  +jgi{m,  n).  (5) 

Let  g\m,n)  be  the  output  pixel  in  the  reduced  sidelobe  image,  and  let 


g(m  -  1,  n)  +  g(m  +  1,  n)  , 

g(m,  n-  1)  +  g(m,  n  +  1)  ,  and  (6) 

g(m  -  1,  n  -  1)  +  g(m  +  1,  n  +  1)  +g(m  -  1,  n  +  1)  +  g(m  +  1,  n  -  1)  . 

At  pixel  (m,  n)  the  sequential  convolution  of  the  1-D  IPR  shown  in  equa¬ 
tion  (2)  with  the  original  unweighted  image  is  given  by 


g(m,n)  =  g(m,n)  +  wnQnandg\m,n)=  g{rn,n)  +  wmQm  . 

At  pixel  (m,  n),  the  convolution  of  the  2-D  IPR  shown  in  equation  (4)  with 
the  original  unweighted  image  is  given  by 

g\m,  n)  =  g(m,  n)  +  wmwnPm>  n  +  wmQm  +  wnQn.  (8) 

We  wish  to  minimize  | g'(m,  n)\2  as  a  function  of  (wm,  wn)  subject  to  the 
constraints  0  <  zvm  <  V2  and  0  <  wn  <  V2.  These  constraints  force  the  window 
to  lie  between  a  uniform  window  (iv  =  0)  and  a  Hanning  window  ( w 
=  0.5).  Without  these  constraints,  the  output  pixel  value  would  always  be 
driven  to  zero.  With  this  background,  we  now  derive  optimizations  for 
four  cases.  We  will  use  I,  Q  in  the  headings  to  mean  the  in-phase  and 
quadrature  components,  which  are  the  real  and  imaginary  components, 
respectively. 


g'(m,  n)  = 


2.1  Case  1:  I,  Q  Separately  and  2-D  Coupled  Weights 


Although  this  derivation  is  provided  by  Stankwitz  and  others  [4],  it  is 
repeated  here  for  completeness.  Consider  the  2-D  case  where  zvm  =  wn  =  w 
(i.ev  the  weightings  in  the  two  dimensions  are  coupled),  and  where 
g(m,  n)  is  either  the  real  or  the  imaginary  component  of  the  image.  The 
operations  are  carried  out  separately  on  the  real  and  imaginary  images. 
Let  Qm  n  =  Qm  +  Q„,  where  Qm  and  Q„  are  defined  in  equation  (6).  In  this 
case,  equation  (8)  becomes 

g'{m,  n )  =  g(m,  n)  +  zv2Pm/  n  +  wQnh  n  .  (9) 


dg\m,n) 

Local  extrema  in  g2(m,  n)  occur  at  the  values  of  w  for  which  — — 


=  0,  which  are  given  by 

—  Qm,  n  —  \j  Qm,  n  ~  ^  m,  n  *0 


2  P 


w  = 


m,  n 


-Q 

2  P, 


m,  n 


;  for  Qm,  n  -  4Pm,  n  g(nl’  n)  >  0 
otherwise  . 


(10) 


Noting  that  g'(m,  n)  =  w2Pm  n  +  wQm  n  =  -Q2I%  „/4Pm  n  for  the  otherwise 
case,  the  output  pixel  is  then  computed  by  the  following  rule: 


1  0 ; 

J  Q  —  Qm,  n.1  ^m,  n  ' 

|  giptf  Zl)  +  ZV  Pm,  n  +  wQmr  n/ 


if  there  are  two  zv  and  either  lies  in  the  interval  [0,  V2 1 , 

if  there  is  a  single  zv  in  the  interval  [0,  V2] , 

where  zv  =  arg  min  ( |  g{m,  n )  +  zvmPnh  n  +  zvQVh  „  |) ,  and 


y=  {0,0.5}; 


otherwise  . 


(11) 


2.2 


Case  2:  l,  Q  Separately  and  2-D  Uncoupled  Weights 


By  allowing  wm  wn,  the  two  dimensions  can  be  uncoupled.  The  deriva¬ 
tion  of  the  solution  follows  the  same  steps.  Again,  the  operations  are 
carried  out  separately  on  the  real  and  imaginary  images.  Local  extrema  in 


g'2(m,  n)  occur  at  the  values  of  zvm  for  which 


d  g'2(m,  n) 


or  g'2(m,  n)  =  0, 


which  is  given  by 


{zvm,zv„)  =  lT^, 

X1  m,n 


Qn  Qm 

P,„  „ 


and 


(™nQn  +  g(™,  n ))  Qr  ^  =  ~(zvmQm  +  g(m,  n )) 


71)  P  +0  n  ZV  P  +0 

m,n  1  >cm  vwm±m,n  1  >cr 


(12) 

(13) 
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Solving  equation  (13)  for  the  possible  endpoints  of  a  curve  passing 
through  the  constraining  box  [0. .  .V2,  0. .  .V2\,  we  get  a  set  J1  defined  in 
equation  (14).  Tossing  points  outside  the  constraining  box  results  in  set  J2- 
If  there  is  a  sloping  valley  running  across  the  constraining  box,  J2  will  not 
be  empty.  If  J2  is  empty,  then  the  corner  points  must  be  tested  to  find  the 
minimum.  We  construct  a  set  /  containing  the  test  points  of  interest  as 


r  l(-(2g(m,n)  +  Qn)  i; 

/l=U  P  +  2Qm 


V  ><m  1\  m 


-(2 g(m,  n)  +  Qm)\\ 

P  +  2Qn  II 


J2  =  J i  n  {O...V2} 


(14) 


r  1  J2r  if  Ji  ^  {empty} 

1  \{(0,  0),  (0,  0.5),  (0.5,  0),  (0.5,  0.5)} ;  otherwise  . 

Now  we  can  solve  for  the  output.  If  from  equation  (11)  g'(m,n),  then  no 
more  work  is  necessary.  Otherwise,  if  (wm,  zvn)  from  equation  (12)  are  both 
in  the  interval  [0,  V2],  then  they  can  be  used  directly  in  equation  (8).  If  they 
are  not,  then  the  set  /  from  equation  (14)  is  computed.  The  output  point  is 
then  computed  as 


g'(m,  n)  =  g(m,  n)  +  zvmzvnPmn  +  zvmQm  +  zvnQn  , 
where  (zvm,  zvn)  =  arg  min  (|  g(m,  n)  +  zvmzvnPm/  n  +  zvmQm  +  zvnQn  |)  . 


(15) 


2.3  Case  3: 1,  Q  Jointly  and  2-D  Coupled  Weights 

Now  we  let  zvm  =  zv  and  compute  the  optimum  weight  for  the  complex 
data,  instead  of  separately  operating  on  the  real  and  imaginary  parts.  We 
will  drop  m  and  n  for  compactness.  Subscripts  r  and  i  indicate  the  real 

and  imaginary  parts,  respectively.  Local  extrema  in  |  g\m,  n)  |2  occur  at 

3  |  g'(m,  n)  |2 


the  values  of  for  which- 
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=  0  ,  or \g\in,  n)\  =  0.  This 


results  in  the  following  cubic  equation: 


3,3  PrQr  +  PiQi  2  2  (Prgr  +  Pigi)  +  (Q2r  +  Ql 

ZV6  +  ^  — — r - —L  ZVZ  +  — - 7 - - - r - ZV 

2  _)_  pr  nt  rt-  ,  n2) 


2  [Pt  +  Pf 


Qrgr  +  Qjg  i 

2  (P2  +  P2 


=  0 


(16) 


This  equation  can  be  put  into  the  form  x3  +  ax  +  b  =  0  by  making  the 
substitution 


zv  =  x  ■ 


PtQr  +  PiQi 

2  [P?  +  P?) 


,  where  a  = 


2 {Prgr  +  Pfgi)  +  ( Qr  +  Q?)  I  PrQr  +  PiQ 


2[P}  +  Pf) 


6(P2  +  P,2 


,  and 


i  )  ) 


b= 


PrQr  +  PiQi) |3  ,  (PrQr  +  PiQi)  (2{Prgr  +  Pigi)  +  (Q?  +  Qi?))  ,  Qrgr  +  Qigj 

2(P2  +  P2) 


P2  +  P-2 
V  1  r  +  1  i  ) 


4(P2  +  P2) 


2^2 


(17) 
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For  the  solution,  we  let 


c  = 


27 


,  and  d  = 


and  the  solutions  for  zv  are 


(18) 


w  = 


c  +  d  — 


PrQr  +  PiQi 


,  PZ 

r  ±  i 


2(P,2 

C  —  d  rjj  _  c  +  d  _  PrQr  +  P,Q  i 
2  2  2 (P2  +  P2) 


and 


c-d 


•/—< 3  ■ 


c  + 


d  PrQr  +  PiQi 


2(Pj 


,  PZ 

r  '  1  i 


(19) 


The  following  conditions  apply  to  the  roots. 

If  (ll  +  fP  >  q  ,  there  will  be  one  real  root  and  two  conjugate  imaginary 
4  27 

roots. 


If  (ll  ,  fP  _  0  ,  there  will  be  three  real  roots,  of  which  at  least  two  are 
4  27 


equal. 

If  P2  +  <  0  /  there  will  be  three  real  unequal  roots. 

4  27 

We  construct  a  set  /  as  follows: 


/i  =  {all  real  roots}  n  (0  . .  .  V2} 


/  = 


//1; 


if  J1  ^  {empty} 


(20) 


\{(0, 0),  (0, 0.5),  (0.5, 0),  (0.5,  0.5)};  otherwise  . 

The  complex  output  point  is  then  computed  with  the  complex  data  using 
g\m,  n)  =  g(m,  n)  +  zv2%,  +  wQ ,  (21) 


where  zv  =  arg  min  (|  g(m,  n)  +  iv2Pm>  n  +  zvQ  |)  . 


2.4  Case  4:  l,  Q  Jointly  and  2-D  Uncoupled  Weights 

Now  we  let  zvm  ^  zvn  and  compute  the  optimum  weight  for  the  complex 
data,  instead  of  separately  operating  on  the  real  and  imaginary  parts. 

Local  extrema  in  |  g\m,  n)  |2  occur  at  the  values  of  zvm  for  which 

d \g\m,n)\  d \g\m,n)\  /  2 

- =; - =  0  and - =; - =  0  ,  or  |  g\m,  n)\~  =  0.  Unlike  the  other 

dzvm  dzvn 
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cases,  the  resulting  simultaneous  equations  could  not  be  solved.  Conse¬ 
quently,  an  iterative  algorithm  is  used  to  find  the  minimum.  The  solution 
to  the  two  partial  derivatives  results  in 


wr 


_  (^n(PfQnr  +  P iQni )  +  ^njPrSr  PjSi  QmrQm  +  QmiQni )  +  QmrSr  QmiSi) 
^vn(Pr  Pi  )  "t  n(PrQmr  PiQmi)  Qmr  Qmi 

and 

_  —  famjPrQmr  PjQmi )  ^vm{PrSr  QmrQm  QmiQni )  QnrSr  QniSi 

wl(P?  +  P?)  +  2  wm{PrQnr  +  P&J  +  Q2r  +  Q2, 


(22) 


(23) 


The  solution  for  |gV/  n)  1 2  =0  can  be  put  into  the  form  aw  2  +  /ray  +  c  =  0, 
where 

a  =  w}n{P 2  +  P2)  +  2 wm{PrQnr  +  PiQJ  +  Q2r  +  Q2-  , 

b  —  2(wm(PrQmr  +  PjQmj)  +  Wm[Prgr  +  Pjgj  +  QmrQnr  QmiQni )  QnrSr  QniSi )  '  (24) 

c  —  ^ m(Qmr  Qmi)  2^ m(QmrS r  QmiSi )  Sr  Si  • 


This  solution — depending  on  the  relative  size  of  a,  b,  and  c — has  the  roots 
of 


Pi  = 


sgn(b)'Jl? 

la 


■4ac 


and  R9  = 


aRi 


or 


(25) 


Ri  = 


-2c 


-b  -  sgn(f>)\/f,2-4flc 


and  R9  =  — 

z  flRi 


(26) 


For  example  if  c  ~  0,  then  equation  (25)  is  better  conditioned.  But  if  a  ~  0, 
then  equation  (26)  is  better  conditioned  and  R2  can  be  ignored.  Areal  root 
is  required  for  g'(m,n)  =  0.  Thus,  for  zvn  in  the  interval  [0  .  .  .  Z2\,  if  b2  -  4 ac 
>  0  (i.e.,  real  roots)  and  either  root  is  in  the  interval  [0  .  .  .  V2\,  then  g'(m,  n ) 
=  0.  Otherwise,  use  equation  (23)  to  find  an  interval 

Hn  =  [ 0  . . .  V2,  raj  n  [0  . . .  V2,  0  . . .  V2\  .  (27) 


Use  equation  (22)  to  find  an  interval 

Hm  =  [wn,  0  . . .  i/2]  n  [0  . . .  V2,  0  . . .  V2\  .  (28) 

If  both  exist,  then  let 


J={(ram,ra„)}  =  HmuHn,  (29) 

where  /  is  now  the  optimum  weights.  If  only  Hn  exists,  then  let 

/  =  {(0,  ra„),  (V2,  raj}  e  H„,  .  (30) 

If  only  Hm  exists,  then  let 

J={(zvm,0),(zvm,V2)\<=  Hm  .  (31) 
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If  neither  set  exists,  then  let 

/=  {(0,0),  (0,0.5),  (0.5,0),  (0.5,  0.5)}  . 
The  solution  is  then  found  with  equation  (15). 


3.  ASR  Algorithm 

The  ASR  algorithm  takes  SVA  to  higher  order  and  uses  a  different  con¬ 
straint.  To  go  to  higher  order,  we  let  the  single  weight  wm  grow  to  Nm 
weights  with  the  notation  kwm,  where  k  =1  .  .  .  Nm.  Similarly,  we  let  the 
single  weight  zvn  grow  to  Nn  weights  with  the  notation  pwn,  where  p  =  1  .  .  . 
Nn.  ASR  adopts  the  constraint 

Nm  Nn 

ZVm+£P^<c,  (33) 

k  =  1  p  =  1 

where  c  is  chosen  arbitrarily  by  the  user.  In  other  words,  if  equation  (33) 
is  met  where  the  minimum  occurs,  then  the  minimum  is  used.  If 
equation  (33)  is  not  met  where  the  minimum  occurs,  then  a  new  mini¬ 
mum  is  found  where  the  sum  of  the  squares  of  the  weights  exactly  equals 
c,  and  this  minimum  is  used.  Like  the  SVA,  if  no  constraint  was  used,  then 
the  weights  would  always  force  the  output  to  zero.  This  new  constraint, 
however,  has  no  good  physical  basis — it  is  just  mathematically 
convenient. 

Along  with  Behboodian,  Reise,  and  McClellan  [9],  we  have  found  that  the 
application  of  the  ASR  algorithm  to  UWB/UWA  SAR  data  is  not  as 
helpful  as  SVA  because  ASR  attacks  large  targets. 

3.1  Physical  Explanation  for  Poor  Performance 

There  is  a  simple  physical  reason,  never  discussed  in  the  literature,  why 
higher  order  SVA  and  ASR  would  have  difficulty  attacking  large  targets. 
Namely,  large  targets  are  not  isotropic.  Instead,  they  typically  have  large 
backscatter  over  a  few  degrees  of  angular  space.  This  is  due  to  the  specu¬ 
lar  response  off  the  side  of  the  target  as  the  SAR  passes.  The  derivation  of 
SVA  assumed  an  isotropic  scatterer,  and  set  up  the  weights  such  that  the 
isotropic  target  could  not  be  attacked.  With  a  first-order  filter,  the  freedom 
to  attack  a  nonisotropic  scatterer  is  very  small.  The  results  bear  this  out. 
With  higher  orders,  however,  and  especially  with  the  constraint  used  with 
ASR,  there  is  great  freedom  to  attack  large  targets.  While  this  paper  does 
not  cover  derivation  and  implementation  of  another  constraint,  we 
recommend  that  future  work  be  directed  along  four  approaches. 

3.2  Proposed  Solution  for  the  Poor  Performance 

We  propose  four  ways  to  mitigate  the  poor  performance.  First,  the  present 
constraint  could  be  weighted,  as  in 
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Here,  W  could  be  made  monotonically  increasing  to  force  the  higher 
order  terms  to  smaller  values.  This  would  be  consistent  with  known  good 
windows  such  as  Nutall  windows. 


Second,  the  present  constraint  could  be  modified  to  limit  how  far  the 
solution  can  move  away  from  a  known  good  window.  For  example,  let 
f  =  [by, 2  wif  3  Wj,  ...  }  be  the  ideal  coefficients  for  a  Nutall  window,  and  let 


W  = 


!jl 


i 


i 


\  107°  '  2Wa  ' 

constraint  would  be 


7 


where  a  good  a  is  probably  2.  Then  the 


E  [kwm  -  kw{)2  Wk  +  E  (pzvn  -  PiVi)2  W  <c  . 

k  =  1  V  =  1 


(35) 


Third,  the  present  constraint  could  be  simply  replaced  in  the  computer 
with  a  nonlinear  iterative  procedure  that  constrained  each  term  to  a  range 
from  0  to  a  limit;  for  example,  (3y/,  where  /3  is  slightly  larger  than  1.  This 
would  be  most  like  the  first-order  SVA  that  constrains  the  solution  to  fall 
between  uniform  and  Hanning  weighting. 

Fourth,  as  will  be  seen  in  the  results  section,  application  of  SVA  to  an 
image  that  was  previously  weighted  effectively  raises  the  order  of  the 
overall  filter.  For  example,  two  three-tap  convolutions  can  be  duplicated 
with  a  single  five-tap  convolution.  Therefore,  SVA  could  be  applied  to 
cases  of  the  original  data,  data  modified  by  a  Hamming  window,  and 
data  modified  by  a  Blackman  window.  Then  the  overall  minimum  can  be 
taken  as  the  output  pixel. 

These  four  approaches  would  limit  the  ability  of  the  variable  apodization 
window  to  place  a  dip  in  the  region  where  a  large  target  has  its  high 
specular  lobe.  At  the  same  time,  they  would  allow  a  higher  order  window 
to  be  applied.  Therefore,  better  sidelobe  performance  would  be  expected. 


3.3  Taking  Advantage  of  the  Problem  for  ATD  Purposes 

The  problem  could  be  turned  around  for  use  in  an  automatic  target 
detector.  In  general,  manmade  targets  have  flat  surfaces  that  produce 
specular  lobes.  Therefore,  the  window  could  be  constructed  specifically  to 
allow  the  flexibility  to  produce  a  dip  where  the  response  from  a  large 
target  could  fall.  For  example,  ()  could  be  changed  to 

A(n)  =  1  +  2 w x  cos  +  #i)  +  2 zv2  cos  +  fl2)  +  •  •  •  +  2wk  cos  n  +  8n  j  .  (36) 

A  detector  could  then  be  trained  (e.g.,  a  neural  net  or  a  polynomial  detec¬ 
tor)  to  observe  the  coefficients  and  generate  detections. 
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4.  Results  on  Actual  SAR  Data 

We  evaluated  the  algorithms  using  data  from  the  U.S.  Army  Research 
Laboratory  (ARL)  BoomSAR  system.  The  image  chip  selected  for  analysis 
contains  a  corner  reflector  as  well  as  vehicles  and  trees,  which  provide 
canonical  IPRs.  Figure  2  shows  the  image  created  using  uniform  weight¬ 
ing.  Figure  3  shows  the  image  created  using  Hamming  weighting.  These 
provide  a  baseline  and  show  the  sidelobe  reduction  and  resolution  reduc¬ 
tion  available  with  a  classic  window.  Figure  4  shows  an  image  after 
applying  SVA  to  the  uniform  windowed  image.  It  shows  that  very  signifi¬ 
cant  improvements  are  gained  through  SVA.  The  resolution  is  equal  to 
figure  2,  yet  its  sidelobes  are  lower  than  figure  3.  Figure  5  shows  an  image 
after  applying  SVA  to  the  Hamming  windowed  image.  As  expected,  the 
resolution  does  not  equal  that  of  figure  3.  Nonetheless,  there  is  noticeable 
improvement  in  the  sidelobes.  The  sidelobe  reduction  provides  greater 
detail  in  the  resulting  image.  Figure  6  shows  the  effect  of  applying  the 
ASR  algorithm  to  the  Hamming-weighted  image.  Figure  7  is  a  plot  show¬ 
ing  the  cross-range  profiles  for  four  cases:  Hamming,  uniform,  SVA-on- 
Hamming,  and  SVA-on-uniform.  Figure  8  is  a  plot  showing  the  down- 
range  profiles  for  the  four  cases.  Figure  7  does  not  fully  capture  the 
effectiveness  of  the  SVA  algorithm  because  the  cross-range  sidelobes  lie  at 
an  angle  to  the  down-range,  cross-range  grid  (i.e.,  they  appear  as  an  "X" 
in  fig.  2).  Nonetheless,  the  figures  show  that  the  SVA  has  reduced  the 
sidelobes  by  over  10  dB,  relative  to  classic  Hamming  weighting. 

It  is  interesting  to  note  that  the  SVA  applied  to  a  Hamming-weighted 
image  sometimes  outperforms  the  SVA  applied  to  a  uniformly  weighted 
image.  At  first,  this  occurrence  may  seem  unexpected  because  the  SVA- 
on-uniform  "should  have  been  optimal."  However,  by  prewindowing 
with  the  Hamming  weight,  the  SVA  is  now  effectively  applying  a  five-tap 
filter  to  the  uniform  image  instead  of  a  three-tap  filter.  This  result  argues 
that  a  higher  order  filter  will  further  improve  performance. 

As  mentioned  previously,  the  ASR  algorithm  allows  a  higher  order  filter, 
but  the  cost  function  for  ASR  appears  to  have  problems.  We  found,  along 
with  Behboudian  and  others  [9],  that  large  targets  were  distorted  by  the 
ASR  algorithm.  Object  size,  however,  is  one  of  the  important  features 
used  in  automatic  target  detection  (ATD)  for  FOPEN  targets.  Application 
of  the  SVA  to  the  Hamming-weighted  image  shows  the  utility  of  the 
fourth  method  proposed  in  section  3.2  for  mitigating  ASR's  attack  on 
large  targets. 
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Figure  2.  Corner 
reflector  image  using 
uniform  weighting. 


Figure  3.  Corner 
reflector  image  using 
Hamming  weighting. 
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Figure  4.  Corner 
reflector  image  after 
applying  SVA  to  the 
uniform  windowed 
image. 


Figure  5.  Corner 
reflector  image  after 
applying  SVA  to  the 
Hamming 
windowed  image. 
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Figure  6.  Corner 
reflector  image  after 
applying  ASR  to  the 
Hamming  windowed 
image. 
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Figure  8.  Down- 
range  profiles  for 
Hamming,  uniform, 
SVA-on-Hamming, 
and  SVA-on-uniform 
windows. 
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5.  Conclusions  and  Future  Work 

We  have  derived  optimizations  for  four  2-D  cases  of  first-order  SVA,  for 
cases  for  real  and  imaginary  parts  both  separately  and  jointly,  and  for 
each  of  these  with  the  2-D  weights  both  coupled  and  uncoupled.  We 
reference  material  that  indicates  that  SVA  needs  no  additional  complica¬ 
tions  to  operate  in  the  UWA  domain,  consistent  with  FOPEN  SAR.  Fur¬ 
thermore,  we  have  shown  that  SVA  is  very  effective  on  UWB/UWA 
images. 

We  obtained  over  10  dB  of  sidelobe  suppression  without  a  reduction  in 
resolution.  The  algorithm  is  capable  of  revealing  to  image  analysts  the 
small  targets  and  signature  aspects  that  were  previously  hidden  by 
sidelobes.  We  expect  the  performance  of  automatic  target  detection  and 
recognition  (ATD/R)  algorithms  to  significantly  improve  after  application 
of  SVA  because  the  analysts  will  have  access  to  information  in  images  that 
was  previously  buried  by  sidelobes.  If  unsmoothed  source  data  are 
available,  then  SVA  can  provide  even  greater  gains  because  previously 
unresolved  features  may  be  resolved.  Improved  detection,  classification, 
and  recognition  by  both  automated  and  human  interpreters  should  result. 

While  ASR  appeared  to  show  promise,  the  cost  function  allows  it  to  attack 
large  targets.  This  action  has  deleterious  effects  on  current  FOPEN  ATD 
algorithms.  Most  targets  of  interest  for  FOPEN  applications  are  large. 

We  produced  a  physical  explanation  for  this  behavior.  Based  on  these 
physics,  we  know  that  the  key  is  in  the  cost  function  constraints.  We 
proposed  four  alternative  constraints  that  would  limit  the  variable 
apodizations  to  be  closer  to  known  high-performance  windows,  such  as  a 
Blackman-Harris  for  second  order  and  Nutall  for  higher  order. 

Unexpectedly,  based  on  the  physical  explanation  for  the  ASR  algorithm 
for  attacking  large  targets,  we  propose  that  advantage  could  be  taken  of 
this  "drawback."  We  propose  using  it  for  ATD  purposes.  In  this  case,  the 
form  of  the  apodization  window  is  specifically  constructed  to  allow  dips 
to  occur.  The  window  coefficients  (as  opposed  to  the  image)  are  then  fed 
to  a  trained  detector  as  one  part  of  the  ATD  processes. 
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it  simultaneously  reduces  the  sidelobes  commensurate  with  the  order  of  the  filter  (i.e.,  the  number  of 
cosine  terms  used).  Results  are  shown  on  imagery  from  the  U.S.  Army  Research  Laboratory's  ultra- 
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